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ABSTRACT 

We find the general solution of equations of motion for self-gravitating 
spherical null dust as a perturbative series in powers of the outgoing mat- 
ter energy-momentum tensor, with the lowest order term being the Vaidya 
solution for the ingoing matter. This is done by representing the null-dust 
model as a 2d dilaton gravity theory, and by using a symmetry of a pure 
2d dilaton gravity to fix the gauge. Quantization of this solution would 
provide an effective metric which includes the back-reaction for a more re- 
alistic black hole evaporation model than the evaporation models studied 
previously. 
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Two-dimensional (2d) dilaton gravity theories have turned out to be very useful 
toy models of black hole formation and evaporation [|T|. They are also relevant for 
four- dimensional (4d) black holes, since a spherically symmetric scalar field collapse 
can be described by a 2d dilaton gravity model 

So = lJ d'x^g [R + 2(V,$)^ + 2e"t>) - ^e-^*(V,/)^ 

where G is the Newton constant and the 4d line element (is4 is related to the 2d line 
element ds by 

(is4' = rfs' + W . (2) 

is a 2d scalar curvature associated with a 2d metric g^y, is the corresponding 
covariant derivative, $ is the dilaton field, / is the scalar field and dVL is a two-sphere 
line element. Quantization of (1) would provide us with a semi-classical metric which 
would include the back-reaction. However, the progress is hindered by the absence 
of explicit solutions of the classical equations of motion. In the case of 2d black 
holes described by the CGHS model 0, the analog of (1) is exactly solvable, and 
by quantizing the solution, one can obtain an effective semi-classical metric up to 
any finite order in matter loops H, ||. Since the matter in the CGHS model is 2d 
conformally coupled, this motivates us to consider a modification of (1) 



^-2 



(3) 



so that / will obey a free-field equation of motion in the conformal gauge, and con- 
sequently the general solution could be found more easily. Action (3) describes the 
dynamics of spherically symmetric self-gravitating null-dust (SSND). For purely in- 
going matter, the solution of the equations of motion is given by the Vaidya metric 
i 

ds^ = -(^l-^^^dv' + 2dvdr , (4) 
where r = exp(— $), G = 1 and 

^-^ = TUv) = l{df/dvr . (5) 

An analytic model of black hole evaporation based on the quantization of the Vaidya 
solution has been studied in [0], and a qualitative agreement with the numerical results 
of has been found. However, in order to fully take into account the back-reaction 
in the analytic approach, one needs the most general solution for which the outgoing 
matter is also present [|, 0]. 
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In order to find the general solution of the SSND model, we rewrite the action 
(3) as 

S=^l£x^[R<P + V{<P)-'^{V,fy] , (6) 

where = e~^* = r^, g^i, = rg^^, and V = 2/r. We do this in order to establish the 
connection with a generic 2d dilaton gravity model, which can be represented by the 
action of the form (6). For example, the CGHS model is given by = 4A^, where A 
is a 2d cosmological constant. The equations of motion are given by 

V^V,(t)-hg^uV = -T^. , (7) 

i?+^ = , □/ = , (8) 
where we have omitted the tildas, □ = g^^V u and 

T^. = d^fd,f~\g,,{y^ff . 

In the conformal gauge ds^ = —e^^dx~^dx~ one obtains 

d+d^<P=~lVe^'' , (9) 

9+V-29+p9+0=-(9+/)2 , (10) 
dJ(l)-2d^pd^(l)= -{d^ff , (11) 

and d^d^f = 0. We want to solve the system (9-12) in analogy to the CGHS case, 
where a gauge p = can be chosen. This is possible since d^d-p = is a consequence 
of the equations of motion. In the SSND case p is not a free field, so that we have 
to find an appropriate modification. We will use as the starting point the free-field 



currents for the / = case [|Tl 



= J 2E + m ' 

j2 = log(V0)2 , (14) 

where 2E = (V0)^ — J{(f)), dJ/dip = V and E' is a constant of motion, j's satisfy 

□ji = , nj^ + R = . (15) 



By going into the conformal gauge one can see that the equations (15) imply that 
ji and j2 = j2 — 2p are free fields. In the SSND case, when / = one has 2ji = 
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r + 2Mlog |r/2M- 1| and j2 = log |r-2M| -2p, where M = -E/4 is the black hole 
mass. The Schwarzschild solution is obtained for 



2ji = U'i^-u) 



J2 = . 



(16) 
(17) 



The gauge choice (16) gives the famihar relation between r and u, v coordinates, while 
the gauge choice (17) is equivalent to 



,2p 



ds^ = dudv = -(1 - 2M/r)dudv. 



:i8) 



Note that the gauge choice j2 = gives a relation between p and 0, which is of the 
type we are looking for, and hence we will concentrate on it. 



When the matter is added, the relations (15) change as 
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+ 



(V0)4 
2 



(V0) 



(19) 



so that j's are not free fields any more. The form of (19) and our remark about j2 
suggest that we look for a free field of the form 



3=J2 + X , 
where X is to be determined from (19). This gives 



d+d-X = - 



(9+0)2 (9_0)2_ 

Equation (21) is vahd for any potential V, and can be solved as 



(20) 



(21) 



X = Xo + loj 



where Xq is a free-field solution. 



d. 



+ 2 dx 



(22) 



In the SSND case cj) = r , = v, x = u, and we choose Xq = 0, which fixes 
the gauge. It is convenient to use a generalized form of the Vaidya metric [10| 



ds^ = -e^^Fdv^ + 2e^drdv , 
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(23) 



which can be related to the conformal form of the metric via ds^ — —C{u,v)dudv 
and C = e^^/r. Then the gauge choice j = is equivalent to 

where Tuu = \T Note that we have not specified the limits of u integration in 

(24), which means that a constant of integration will occur. Only when this constant 
is specified, the gauge will be completely fixed. 

The equations which determine r are (9) and (10), while (11) and (12) are the con- 
sistency conditions for the gauge choice, which are satisfied by construction. Equation 
(9) becomes 

dudy = el'dur , (25) 

while (10) gives 

d,''r^-2d,pdy^-{djf , (26) 
where p is determined from the gauge choice as 

2p = log I - 9„r ■ (27) 

Equation (11) follows from (27) and (24), while the equation (12) follows from 
(27), (25) and (24). By integrating (25) with respect to m, we obtain 



(28) 



m{v) is a constant of integration, and can be determined from (26). By inserting (28) 
into (26), and by using (27), we get 

dTfi 1 1 r 

" ~ 2 ~ '''^ ~ 2r9,V')9,r - -d, j du{e^ - l)d^r . (29) 

When ■0 = 0, (28) yields F = 1 — 2m/ r and (29) gives the relation (5), so that 
one recovers the Vaidya solution. When 7^ 0, the equations look difficult. However, 
their form is such that a perturbative solution in T^^ can be easily found. Let us 
introduce an expansion parameter e by replacing r„„ with eT^u- Then we will seek a 
solution in the form 

r = ro + eri + eVs H . (30) 

The expansion (30) then implies 



-0 = e'^i + e^V2 H 
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(31) 



where 

-01 = / du "5'" , -02 = - / duTuuP^^^^^-^ , ■■■ ■ (32) 
In general one should also take 

m{v) — mo{v) + e'mi{v) + e^m,2{v) + ■ ■ • , (33) 

although it is possible that the series in (33) gets truncated, like in the example we 
consider in this paper, where m = euii. By inserting the expansions (30), (31) and 
(33) into the equation for r (28), we get an infinite hierarchy of equations 

OvTi = -^n ^ 7^ / duipidurQ , 

ro^ ro 2ro J 

dvr2 = ^r2- — -—dyri + f dulU'ipi'^ + 'ip2) duTQ + 'ipiduri] , .... (34) 

The system (34) can be solved, since for every n the equation for r„ does not involve 
rjt with k > n and each equation is a first order linear differential equation for r„, 
except for n — 0, which is a non-linear first order differential equation. Therefore by 
starting from n — 1, one can write an explicit solution for any r„ in terms of ro and 
Tuu- At each step several integration constant C{u) and C{v) will arise (due to u and 
V integrations), and these constants can be determined from the constraint equation 
(29), exactly as in the T„„ = case. This can be done because the equation (29) also 
decomposes into an infinite hierarchy of equations under the expansions (30), (31) 
and (33) 

T„„ = drrio/dv , 
= dmi/dv + + 2ro9„V'i)(?^ro - ^dy j dui/jiduVo , 
= dm2/dv + ^{ijj2 + + 2rodyijj2 + 2ridyijji)dyro + ^(V'l + 2rodvijji)dyri 
- ^dy j du[{'4j2 + h'4'i'^)duro + i'lduri] , ... . (35) 

The equations (35) will also determine the m„(v), provided we fix the integration 
constants ipniv) which appear in (32). These are related to the complete specification 
of the gauge, or equivalently, to the choice of the v coordinate, since a coordinate 
change v — v{v) in (23) gives 

■0 = + log ^ . (36) 
av 
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At the end, one sets e — 1 and writes the solution as 

r = ro + ri + r2 H , V = V'l + V'2 H ■ (37) 

These general features can be nicely illustrated on the example of a white hole 
emiting a shock-wave. This is simply a shock-wave Vaidya solution where the u and 
V coordinates are interchanged. The corresponding spacetime is described by the 
Penrose diagram of Fig.l. In this case ^ 0, and one can find an exact solution 



Figure 1: Penrose diagram of a white hole emiting a shock- wave. 



for r = r{u, v), so that the expansions (30) and (31) and the corresponding equations 
(34) and (35) can be checked. By taking T^u — M5{u — uq), we get 

—du^ — 2dudr u > uq 
-(1 - 2M/r)dv? - 2dudr u < Uq 



ds' = 

This can be rewritten as 
ds^- 

where F = 1 - 2M/r, 



—dv^ + 2dvdr u > uq 
-Fe'^^dv'^ + 2e^dvdr u < Uq 



1 - 



4M 



V -Uo 



(38) 



V' = - log 

and r — ^{v — u) ior u > Uq, while for u < Uq 

r + Tglog \r/rs — 1| = \{y — u) + r^log \{y — Uo)/2rs — 1| , 
where Vg — 2M. The form of the solution, given by (40) and (41), is such that 

n>0 n>l 



(39) 



(40) 



(41) 



(42) 



which is of the form (37). More exactly, one can show that r„ = r^^f^ and ipn = Vs^il^n 
satisfy the equations (34) and (35) with m = mi = M. Also, starting from n = 2, 
non-trivial integration constants ipniv) appear. 

In conclusion we can say that we have found a useful form of the general solution, 
given by the expansions (37) and the equations (34) and (35). The expansions (37) 
are clearly in powers of the outgoing energy-momentum tensor, and by truncating 
them at finite n we obtain an explicit perturbative solution. This form of the solution 
can be used to construct an effective quantum metric in the approximation of a 
finite number of matter loops via the method of quantization of the classical solution 
0, ^, 1^, 0] . The corresponding construction is going to be more involved than in the 
2d case ^, Q or the Vaidya case |0, since the relation between the metric and the 
dilaton (radius) is more complicated. For example, for the one-loop approximation 
one would take r = tq + ri and ip = tpi, with T^^ replaced by (T^;/) evaluated in 
an appropriate quantum state. Note that in the one-loop case one does not have to 
truncate the expansions (37) at n = 1. By including the higher-order terms, one 
extends the validity of the one-loop approximation to a smaller radius. 

Our solution can also serve as a good starting point for finding approximate 
analytic solutions for the more realistic collapse described by the action (1). 
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